{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Housekeeping"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "code_folding": [],
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Bunch of imports \n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "%matplotlib inline\n",
    "import pandas as pd"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "colorstyle=dict(red_face = np.array([255,85,65])/255, red_edge = np.array([201,67,52])/255,\n",
    "                blue_face = np.array([49,115,255])/255, blue_edge = np.array([36,85,189])/255,\n",
    "                green_face= np.array([84,224,81])/255, green_edge= np.array([62,166,60])/255,\n",
    "                yellow_face=np.array([255,207,49])/255,yellow_edge=np.array([191,155,36])/255,\n",
    "                gray_face=np.array([169,169,169])/255,gray_edge=np.array([137,137,137])/255)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {
    "collapsed": true
   },
   "source": [
    "# Load saved profile"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import scipy.io as sio"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [],
   "source": [
    "Data = sio.loadmat('FigS2B.mat',squeeze_me=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "TTilde_2kHz=Data['TTilde_2kHz']\n",
    "TTilde_2kHz_err=Data['TTilde_2kHz_err']\n",
    "TransRate_2kHz=Data['TransRate_2kHz']\n",
    "TransRate_2kHz_err=Data['TransRate_2kHz_err']\n",
    "\n",
    "TTilde_5kHz=Data['TTilde_5kHz']\n",
    "TTilde_5kHz_err=Data['TTilde_5kHz_err']\n",
    "TransRate_5kHz=Data['TransRate_5kHz']\n",
    "TransRate_5kHz_err=Data['TransRate_5kHz_err']\n",
    "\n",
    "TTilde_fit=Data['TTilde_fit']\n",
    "TransRate_5kHz_fit=Data['TransRate_5kHz_fit']\n",
    "TransRate_2kHz_fit=Data['TransRate_2kHz_fit']"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Fig2B"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAJsAAACICAYAAADwI5m5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAFbNJREFUeJztnXl41NW5xz8nyWSfhJCEbOxl3yEB\nxdomoqCo2FVrUVuXir2IWgW5t7d6i/b23j4u97Z1p27VikqvthUasHUBVxAJBMJiZAuQZLIvk8ks\nycx7/5gkBsgkM8nMbyaZ3+d55nkyOed3zvkl3+c92/ueo0QEHR0tiAh2A3TCB11sOpqhi01HM3Sx\n6WiGLjYdzdDFpqMZuth0NEMXm45mRAW7Ab2RlpYmY8eODXYz/IrZbMZoNAa7GQFh9+7dtSKS7ik9\npMU2duxYPv/882A3w68UFBSwbdu2YDcjICilynpL17tRjRmqQvMGXWwas2LFimA3IWjoYtOYvLw8\nTevbcqxQ0/p6Qxebxmht2bYe36Jpfb2hi01jtJyJmiwmGmwNmCwmzersDV1sGlNRUaFJPSaLiVXb\nV1I3vI5V21eGhOB0sWlMIGajZ4/LthwrZO22NVhTrMRnx2FNsbJ225qgj990sWnM+vXr/VqeyWLi\n1UMbzrBcS8dfzkMFjxDXEEdrhZW4hjgeKniEpeMv92vdvhLSi7pDkU2bNvmU3+wwU1y9l7e+/Bv3\nLFjDafNpPjz9Aa1trTTZGyk1l2JPtnN94XKmpUwjNiqOS8ZcwpJxl5I7PI8PT3/A3Jx5bD+1jR9M\nuZYi025Ot5zG4XRQbanmyglXkhaXTpO9icToRBIMCUSor2yQyWKipKaEGekzyEzIHNC762LTmOXL\nl7Nhw4Ye07YcK+Si0Yuoaa1mVNJo1u99mj3VezjWfAxXootV21fyy/kPsHjMYoqqiiip3Y8j1UFC\ndjytEVYyI7O4bc5PiY2KBeCG6TeQlZBFXmYeUREGAKKjYmi2N7Oh9BXaYtt4t/Id7pnj7mJbHC1Y\n263cNPMm5mcu4Nc7fsUnlZ/gSnASWxzH+kuexdJmIS4qjuzEbJ/fXRebxoxcMLLH3++s2MlTe57k\nb0f+ytdzLuSG6T8iPX4ELY4WonKiiM92d4n/OPE2q+evYUb6TJaMu5RV21diVVbim+O5Of8WUuNS\nv6rLOIobZ950Rj2nmk/yjxNvozIVw7KTaa2w8szep0iJTeGycUtZMvZSBGHr8S0cqD2ApAgxyTHY\nm+ys3baGiSkTmT1iDldN+JbP766LTWPsM61dP4sItdZanOLkgZ2/pD2pnWpHNYvHLgHgO5O+y8Kc\nC1i1fSWtuMde1+Vf3/V8ZkImj+c/ydpta3io4BGvurml4y9nbsY8Vm1fiU3ZusZzZz97xdeuJDcz\nzy3mOCtGm9HrOjyhQjmULy8vT4bSRrzJYmLW4pns++d+XOLi9cOvYXaYiY+M433X+13WK7U+lR9O\nXd41oDdZTLxy8E9cN+36Hv/Zd717B7+7+DGf23KgtoTpab2PxfqquztKqd0i4nGLRBebRnSue1mM\nFtqr2hlnHMc1k6+lYHQBtdZatwVJcVuvx/Of9MmCbDlWGPSZJvQttoAtfSilblRKfaKU+lgpNc9D\nngeUUkcC1YZQoqSmBEtsC/Z6O4asKMrt5cxIn0FURFRXd7go6mKfhQaEhNC8ISBiU0qlAHcCBcD1\nwO97yJMBTApE/aGAyWLinRPvYLKYEBF2VHyKtdpGdLKBhJEJkM4ZC62ZCZmsnr9mwMsLoUygJggL\ngA9FxAEcV0oZlVIxImLvlud+4L+BNwPUhqDR2WXah9uIORjLoqyLEYQHv/4rHvz8AVpjrB4H5kOZ\nQHWjqUBDt++NwPDOL0qpiUCiiOzrrZCamhry8vLIy8vz+8p7IHnlwMtYU6xEp0VjS7EyPHY4/3be\nz/nmqHzML1r63V0OdgJl2eqBYd2+J3f8rpN1wH/0VUh6evqgdAu/bvoNfPDuBzRWNWGMNFIw6yIM\nke5F1ft+ch/L5i8LcguDQ6DEthP4T6WUAcgCWs7qQscDTyilALKUUr8XkTsD1BbN+bKhlNHxo1mQ\ndR6Lxy45w4IVFBQEr2FBJiBiE5EGpdSTwHZAgLuUUnOAxSLysIgs7MyrlDoylIQmIpSby7lv4X+Q\nkZBxTnp2djZmszkILQs++jqbn2hztvHnLzayeOwS0uM9RrMNaYK2zhZOWNut/L7odzTZG0mOSe41\n72Ca6PgbXWz94Ow1tCeKHifBkMCEYZOot9X3+uxgsdSBQO9GfaT7GpqhNponCp6iwd7A/Tt+gSPV\nTkx9bFgua4DejfqdzjW0qBQDje2NPLX3SSrMFThS7cRmxWJNsfLKwT95fD6cZ6O62Hzkuuk3YKg1\n0HSwiUSVyMq5tzMjfQYx9bHYTQ63G9C06z0+v27dOu0aG2Lo/mw+kpmQyeWjr8AQYWDp+Mu7usvH\n8590u+zM6N1lZ9KkIbsd3Ce62HxgX80+0uLSWDH7tnPSMhMyvRqn5eXlaRbOF2ro3aiX7KrcxR9L\nXsThtPeduRfCVWigi80rdlR8ysbDr3F33t2MTR43oLIeffRRP7Vq8KF3ox7oDGGbnjadjIRMVi+4\n1y/LGeFs2XSx9UDnWlpLtJnIvVE8t/gFv62bhbNl07vRHnjlwMs0RzXjdLpwpbl6XTfzldzcXL+V\nNdjQxdYDczLmIU2CIcFAQnNCr+tmvhLOe6MB60aVUjcCK3C7GN0hIkXd0q4FVgEuoBlYLiLNgWqL\nt4gITY4mLhp9EeOSx/GXL9/kuvy+Q9h8Yage3uwVIuLTB7cz5J3AHUCWhzwpQBEQDYwDPjorPbrb\nzw8Ct/dUTm5urmiFy+WSjYdel//d9T8BrWfixIkBLT+YAJ9LL9rpTze6GtgCvA3c4yFPV8CLiBwH\njEqpmG4Cd3TLmwAc6Ec7/IaIsPHw63zR8AW3zg7syZClpaUBLT+U6Y/YHMCJjo/DQ55eA14AlFK3\nKKX2A9/Ag9i0Cnips9ZxtOkoBhWF2RFYL1p9b9RLlFLJwOvAb3CPxV7ykLWvgBdE5DngOaXUWuBe\nYO3ZhQQ64MUlLnZV7mJU0iiK6ndjTbGya/uusHURCjRei00ptRX3WCwCt3CKxHMoXq8BL0qpWBGx\ndXxtBOL70/iB4BIXLx94iSqLCYSuUxpbcZ/S2P2sDX+iWzbvMIrIeQBKqeGAxwUj6SPgBbhXKXVx\nR/Z64OZ+tb6fuMTFSyV/pMZaw525P6PR3shn2z/r9VQffzFp0qSwHbd57amrlFoHPCsipwPaom4E\nylO3oqWCVw+9ytWTv8/opDGA96f6DJTS0tJB6WZkb62hpe4ALXWHiYpJImfq8nPy9OWp64tl+xaw\nsqM7/Qz3NHeHz60OIi5x8ey+P7DlVCHONCeffvRx1/jMWxehgRIqYXw2czlN1UUkj5hHrDHnjDQR\nF9amE7TUHcLeWk1r0zGiDEZASEydijF9Vr/q9CkGQSkVC8wE5gJzRGRlv2r1En9aNqfLyX0f/jv7\navYhWULCyPgez0ILNLm5uezevVuTujxhM5dzYOs1ZCSbqWoyMvWSl2lvt6AiIjGmTuPAu3ch4iQ2\nMYe642+SM9xGVXMS0y/beI4wuzPg89mUUu8B+4G3gG0i4uznO/qMP8X2pwMvU91azdWTr+Huj+7C\nkeogui56SM88e7JelaVvULbnabKMlYzKNLDncCsWezSR0YmMGH85Exfej7PdTvWxzZws/gM5SeWM\nzDBwuqqN8uYcRs++laxJ3+uxPr8cBqiUmgFcBeQDtcBmoFBEmnz/E3iPP8TmdDlRSmGymEiLSyM6\nMlqz8VlPrF69WhPPj07rNSKpmaomI5Mueob6U9toqTuEw1pDW8txclLtVDTEM+Oy14gfdq6f3tkW\ncKCWzasxm4iUACXAf3Wcq7YMeF4p9TcR8bTWFnScLifP7X+WCcMmsGjMxV2/12p81hPZ2b6fsu0r\np0pe5GTxH0iIqqOuEYbFNXHw/dUMz7mA8Qv+lbikUdhbKikrfoZZF9zmUUCxxhymX7aRpuo9TB8x\nt1eheYPPcaNKqZuAvSKyZ0A1e8FALFun0FrbWrl97qquU4SGGpaGo9QcL0SpSGyWCsbnrcFuMVF1\ndBO1R14nY5iNOktyn1bJH/hzNtrJBOAKpdRLwA4Rqe536wLIx+UfhaTQsrOzffbWtZnLKSt+mjGz\nf0qkIYGWuoO0Np8gdWQ++wuvJjqiFbszjglf/w2RhngSU6eQmDqFnKnLKSt+hunf8Gy9tMTX2Wiy\niDQppc7D7dFxfscibUDoj2VzupzUWWtJi0/H6XKGlNDA7RbubVcqIphr9nH4vVvJTrFysjqS2OSJ\nGNNn42yz0GgqYmRyRdcA3pLwHSZf+GCA38AzfrNsndtVyn2oWj2wR0R+7oc2DpjOeIFpaVPZfHQz\nCsWts1cQERl6vqGlpaU9iq1z5mhMnUH8sHGc2v88ptI3sZpPkpLQxsiMOFwuB62JU5i48L6uZw5s\nvQaU2b00ccG5IYahREC2q7SkM17AlmKlvcjJN7PzuXf+OXv6IcO6devOuJlPRKg/tZ3SD+4mSllp\nc8Uw79uFJGfkkj7uMsTl4uDbP+B0tZmalmSmX/jTrmf9PYAPNINuu2rLscIzbgdusDVQN7yOmNRo\nmo+YyXJlMTzuK2+my8YtDZmj27uvexni0nC0VhFrHMXnby4jnqOMy4mhrrGdCvPIM9az3M/tITnE\nBeW3SzeUUnuAHECz7SpvxmwVLeXc9u4KZISL2H5cWKEVnV1eUkwD1Y1RxCVPJn38UkbNvLkrLTPZ\ngqkpQZOZYyDw25hNROaetV31IyCoe6MucVF4rJCFIxaSlzmfWbNmh6TQAJqqixie0ERlTTuxBidt\nLUdIH3spMPi6w/7i09JHhw/aro5Pr/QR8LIW+B7QjjtW4U7xccHPJS5e2P88ZoeZNQvWEh0Z7cvj\nAefMLnM4Nce30ljdxqiMKEZnxlBe3U7x1lu6ustYY86QFVknAYmu6nbDy/m4u96XgQu7ZfmLiDzU\nkXcjsAh419d6RieNIX9UfkgKrXObp/izGGKSppCcmUf25Gsp+vvNREZGUdWcxOzLnhvyAutOUG54\nEZEvu+W147ZwXuESF68depWC0RexeOxi/7baT5QVP01GUjPpwyKprGlmWPb5jJ17OwBzvr0Zp+XL\nId1deiIoN7x0opTKx+02/kFPhZwd8OISFy+WvEClpfKMS1xDCREhKX02p2oN1DS6UIZUMid8uyv9\n012lZHztyrATGgTvhheUUrNwB84s8zReOzvg5bVDr1JvreeO3DuJiYzp6ZGgYrdUUVb8NO32RqYU\nPI7DVs+MsyzY+vXrWbZMv+HFn/QV8DIBeB74nojUelvoN0flkxqbGnJCExHsLRWc2PMY8ckTGHn+\nL4iI6PlPu2nTJo1bFzoEpBsVkQagM+DlVeBnSqk5Sql7O7L8Frfl+6NSaptS6gpvys1OzCYmKrSE\nZm0+ycH3fsb+v3+XRFshNaUv4LBUecy/fPm5vvvhQsDO+hCR53Fbr+7s7Ui7MlD1aoHNXE6jaReO\n1mrqyz8mIWUCSSmtjBxhAGmmrPgZjxviV145qF99QOjns/mIzVxOyZaryUiycLpGMePSl4iOS+XA\n1q0oLzbEdcum4xUup52D799DpKueUZnxKNVGZen/MfnCB73eAVBK4eP69ZBBF5uXtNQd4kTRYySm\nTqHZepLyausZVszbHYBwFRroYusTZ1srEVGxONttjJzxY4ZlnYfNvKLf+5gbNmwI265UF5sHbOZy\nKkvfoMm0i3Hz7yE5Y25X2kD2MTdv3qyLTecrWpvK2PvWVRgibEhEMjFxI/xW9oYNG/xW1mAj9Pym\ng4TNXI7pyFs0V+3l1P5nSTW2MW9qPNnDrZQVP+O3esJ19wB0ywa4hba/8PsYaKTFFsWMJc9ztPI9\nKmv979u/YkVgT7YMZXSxAV/u+DUuRz3pGQZSXcIXH91P1uSbMCRk+t07I5yvgAxrsdktJlSEgbFz\nVnL43d0IVqrNxoD6mWVnZ4fMSUZaE5Ziszaf4lTJ87TUHmRc3t2kZJ/PzCve1MQtO1yFBmE4QbA2\nn6bor1fQUvYqTmsZccZRgHs5Qws/s3C+dCNsxOZytXNk528o3noL2alO5k6JI2uYheKtt1BZ+oZm\n7Qi1O++1JGBiU0rdqJT6RCn1sVJq3llpFyil9iulbEqpkYFqQyeWhiMc3n4v4nQwa8lT1LcOo7zG\n2RUH4Om8sUAQzpbN5xtevPnQ9w0vyUAisA0Y6akcf9zwYq49KHsLfyy1J98Xl8slIiLW5tNiOrJJ\nrM2nB1y+r+Tn52tep1bQxw0vwQp4aQK3B0SgMNcdRJwOjOmzmHbRbzHEfuWlHsywuXA+mj6oAS99\n4e0NLzZzOVVHN2Ezl+Nst3KyeD3Hdz2KuJwoFXGG0ILNYDwp3F8ELeDFG7y54aV7jOaBIiOJWYuI\niR/BtEW/JSo69G7Ay8vLC9vblIMS8DIQKkvfwFT6Ztd3h7WOrMQmxAXpiU1UVOwkJj6dhgr3yRCZ\nk76r6QSgL8JVaEBgJgjusSI3A58AHwN5wBzg3o60ScA7uLvaD4F/6akMbyYIpi/fkg9emCb7/jxF\ndr2+MCiDfl945JFHgt2EgEGQJgh9BbyUApcMtA67xURt2T+ZtugJ2ttaQv5IKQhvy+bzAc5a0tOR\nWSJC/alt2CyV5ExdjogLpcJmbTqk6evIrEH1X7K31nDk019RdfQthmWdBzDohJabGxIHdgaFQbUR\nX1f2DompU8mY+B2PEeehTjjvIAyq/1j21B8GuwkDxmgMveUYrRhcfVA/CSVr4q+I+FB6J2/RxaYx\n/rrYNpTeyVvCQmyhRDjvjYb00odSqgYo80NRabhvExxKhOI7jRGRdE+JIS02naGF3o3qaMaQElso\neQf7iz7eaa1SamdH2mMqkA6C/qC3jdPB9MFP3sGh9PHinSZ2+3kjcHGw29zbZyhZti7vYBE5DhiV\nUl1noopIk4i0BK95/aKvd+r3Ef/BYCiJzS/ewSGGX474DxUG1XZVH/jFOzjE8MsR/6HCULJsO4EL\nlVIGpdRo/OgdHER6faduR/xfKz4c8R8shtQ6m1LqZuAnuC9nuwv3GGaxiDyslJqE+7j8XKAE2CAi\nTwWtsV7SxzttBqYAnXfAPiwifw9OS/tmSIlNJ7QZSt2oToiji01HM3Sx6WiGLjYdzRhK62whgVLq\nNuCHuF2AhgFHACdwKfBP4DVP6SIS0jsAA0WfjQYIpdRdQLSIPNzxfSawQkTu6Ck9HNC70cAxC9jX\n7ftVwFu9pA95dLEFjrPFlI/b48RT+pBHH7MFAOWOnB4hIpUd3zOBOhFp6ym943c3ArNxu8G/JiIm\nzRseYHTLFhgm4h74d7IM2NxL+tnUBKJRwUa3bIHh7C5yKXBLL+mdvC4iOwLZsGCiz0Y1QCl1o4i8\n2Fce4LAuNh0dP6CP2XQ0QxebjmboYtPRDF1sOpqhi01HM3Sx6WiGLjYdzdDFpqMZuth0NOP/AVRP\nNWZqbjOzAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 144x144 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig = plt.figure(figsize=[2.,2.])\n",
    "ax=fig.add_axes([0.23,0.15,0.75,0.65])\n",
    "\n",
    "plt.errorbar(TTilde_2kHz,TransRate_2kHz,xerr=TTilde_2kHz_err,yerr=TransRate_2kHz_err,marker='o',linestyle='None',\n",
    "            mfc=colorstyle['yellow_face'],mec=colorstyle['yellow_edge'],\n",
    "            ecolor=colorstyle['yellow_edge'],markersize=3,elinewidth=0.75,mew=0.75,label='2 kHz')\n",
    "\n",
    "\n",
    "plt.errorbar(TTilde_5kHz,TransRate_5kHz,xerr=TTilde_5kHz_err,yerr=TransRate_5kHz_err,marker='o',linestyle='None',\n",
    "            mfc=colorstyle['green_face'],mec=colorstyle['green_edge'],\n",
    "            ecolor=colorstyle['green_edge'],markersize=3,elinewidth=0.75,mew=0.75,label='5 kHz')\n",
    "\n",
    "plt.plot(TTilde_fit,TransRate_5kHz_fit,linestyle='--',color=colorstyle['green_edge'],alpha=0.8,linewidth=1,zorder=5)\n",
    "\n",
    "plt.plot(TTilde_fit,TransRate_2kHz_fit,linestyle='--',color=colorstyle['yellow_edge'],alpha=0.8,linewidth=1,zorder=5)\n",
    "\n",
    "\n",
    "plt.axvline(0.167,color='k',linestyle=':',linewidth=1)\n",
    "\n",
    "plt.xlabel(r'$T/T_\\mathrm{F}$',fontsize=9,labelpad=1)\n",
    "plt.ylabel(r'$n_f/n_0$',fontsize=9,labelpad=1)\n",
    "# plt.xlim([0.075,0.175])\n",
    "# plt.ylim([0.0,0.48])\n",
    "ax.tick_params(axis='both',direction='in',labelsize=9,length=3)\n",
    "\n",
    "fig.savefig('FigS2B.pdf',dpi=300)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python [default]",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.3"
  },
  "toc": {
   "colors": {
    "hover_highlight": "#DAA520",
    "navigate_num": "#000000",
    "navigate_text": "#333333",
    "running_highlight": "#FF0000",
    "selected_highlight": "#FFD700",
    "sidebar_border": "#EEEEEE",
    "wrapper_background": "#FFFFFF"
   },
   "moveMenuLeft": true,
   "nav_menu": {
    "height": "30px",
    "width": "252px"
   },
   "navigate_menu": true,
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "threshold": 4,
   "toc_cell": false,
   "toc_position": {
    "height": "630px",
    "left": "0px",
    "right": "1921.3399658203125px",
    "top": "106px",
    "width": "212px"
   },
   "toc_section_display": "block",
   "toc_window_display": true,
   "widenNotebook": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
